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This Commentary on the paper, "Total Hadronic Cross Section Data and the Froissart-Martin 
Bound ", by Fagundes, Menon and Silva (Braz. J. Phys., Vol. 42 (2012); arXiv:1112.4704) was 
invited by the Editors of the Brazilian Journal of Physics to appear directly after the above authors 
printed version, in the same journal issue. We here challenge that paper's conclusions that the 
Froissart bound was violated. We will show that this conclusion follows from a statistical methodol- 
ogy that we question, and will present compelling supplementary evidence that the latest ultra-high 
energy experimental pp cross section data are consistent with a In'^ s behavior that satisfies the 
Froissart bound. 
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Introduction. — We have been invited by the Editors 
of the Brazilian Journal of Physics to write a back-to- 
back Commentary, to be pubhshed in the same printed 
journal, on "Total Hadronic Cross Section Data and the 
Froissart-Martin Bound ", by D. A. Fagundes, M. J. 
Menon and P. V. R. G. Silva, Braz. J. Phys., Vol. 42 
(2012);'arXiv:1112.4704, hereafter referred to as FMS 1]. 
We virish to challenge their conclusion that the original 
Froissart bound [2| is violated, i.e., their statement that 
the total pp (pp) section rises more rapidly with energy 
than cln^ s, where c is a constant and s is the square of 
the hadron-hadron center of mass system (cms). 

In essence, FMS conclude that when all total cross 
section data for pp and pp between 5 < -^/i < 1800 GeV 
are used — a total of 163 datum points — the measured 
total cross sections satisfactorily satisfy the relations 

(JLE = "1 (^7"^ ^"^^(t") ' 

(THE = a + /31n^— , 7 = 2, (2) 

Sh 

CTtot = 0-LE+(yHE, (3) 

where ai, a2, &i, 62, a, 7 and Sh are real constants and 
Si = 1 GeV^ and the -I- sign is for pp and the - sign is for 
pp collisions, i.e., they find that the Froissart bound is 
satisfied, but that is is violated (7 > 2) when one addi- 
tional point at 7 TeV — the Totem Q total cross section 
measurement of pp, titot = 98.3 ± 0.2(stat.) ± 2.8(syst.) 
mb — is considered. As examples, see the 7 = 2 result 
from Table 1, Direct Fit of FMS, and 7 = 2.104 ± 0.027 
from Table 2, VI of FMS. We will show that the statisti- 
cal probability of the 7 = 2 solution with the inclusion of 
the Totem result (Direct Fit model of Table 2, FMS) is 
higher than the statistical probability of the fit (Table 2, 
VI model) that gave 7 > 2. Thus, we find no statistical 
evidence for their claim. 

Further, since they only include one additional high 
energy point at 7 TeV in their analysis, an alternate (and 



in this case, perhaps more transparent) analysis method 
would have been to use the 163 pp and pp cross section 
datum points in the energy interval 5 < ^/s < 1800 GeV 
to predict the 7 TeV total cross section. We will show that 
when we reanalyze the FMS results, the low energy fit 
for 7 = 2 after including the errors of prediction (due to 
statistical errors in the parameters) satisfactorily predicts 
the 7 TeV Totem total cross section measurement. 

We will also comment on the necessity for the simul- 
taneous inclusion of the p data, the ratio of the real to 
the imaginary portion of the total cross section, together 
with the cross section data, i,e., a global fit of p and crtot, 
a statement alluded to (but not carried out) in the Ap- 
pendix of FMS. Since all but one of the parameters is 
common to both p and cr, the simultaneous inclusion of 
the many p- values is important to an accurate fit, in order 
to minimize the (correlated) errors on the fitted parame- 
ters. We strongly disagree with the FMS comments 1 to 
6 in their Section 2.2 concerning the p- values, in which 
they dismiss them as essentially useless, whereas in their 
Appendix, they give a rather cumbersome evaluation us- 
ing their Variant 3 model, to separately evaluate p, and 
not globally evaluate it with the total cross sections, as 
demanded by analyticity. 

Recently, eight measurements of pp cross sections have 
been made at energies higher than the Tevatron energy of 
1800 GeV. At the LHC cms energy of 2.76 TeV, ALICE 
^ has measured cri„oi, at the LHC cms energy of 7 TeV, 
ALICE ATLAS [5|], CMS [gj] and TOTEM have 
measured fTmci; ftot has been measured by TOTEM 0]. 
Most recently, the Pierre Auger Observatory has pub- 
lished pp cross sections for iTinei and crtot at 57 TeV 0, 
using cosmic ray measurements of extended air showers 
to measure crp_air- Most Hkely, this is effectively the high- 
est energy reach that one will ever have experimentally. 

Finally, using these ultra-high energy measurements, 
we will summarize a very recent analysis of Block and 
Halzen [ij that shows that the proton asymptotically be- 
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comes a black disk, whose total cross section asymptoti- 
cally varies with energy as a saturated Froissart bound, 
i.e., as In^ s, and whose asymptotic ratio of inelastic to 
total cross section is 1/2, that of the black disk. This 
result, using all available high energy data, including in- 
elastic cross sections, clearly contradicts the conclusions 
of FMS. 

Statistical Probabilities — For a fit, the goodness- 
of-fit criterion for a given model is normally taken as 
^(Xoi^)j tfi^ integral of the probability distribution in 
for v degrees of freedom, integrated from the ob- 
served minimum Xo to infinity; it is P(xoj ^) that allows 
us to statistically distinguish between models. From Ta- 
ble 1 of FMS, we concentrate on the Direct Fit (DFi) 
and VI models, which give Xo = 145.236, v = 156, 
yielding PdFi — 0.721, whereas for VI, we have Xo = 
145.235, V = 155, yielding Pyi = 0.701. Thus, even 
though FMS relaxed the input restriction that ^ — 2 and 
allowed it to be fit by the data, we get the somewhat 
strange result that FMS have a better, somewhat more 
reliable fit when they fix the value of 7 at 2, the Froissart 
bound limit, than when they allow it to float, suggesting 
perhaps that the true minimum x^ was not achieved in 
their minimization process. In any event, FMS concluded 
that the value 7 = 2 was correct for the energy interval 
5 < < 1800 GeV. 

Next, we calculate the goodness-of-fit probabilities for 
the Direct Fit and VI models from Table 2, which now 
includes one additional point, the Totem point at 7 TeV. 
For the Direct Fit model (DF2), we find xl = 146.01, v ^ 
157, giving Pdf2 = 0.725; this is effectively identical to 
the comparable value obtained (PdFi =0.721) not using 
the Totem point. Thus, if 7 = 2 is satisfactory for the low 
energy data, it appears to be exactly the same level of 
confidence when we include the Totem point. Finally, we 
calculate from Table 2 the high energy version of the VI 
model (VI2), i.e., Xo = 145.86, v = 156, giving Pyia = 
0.709, a probability smaller than that of the Direct Fit 
model (DF2), i.e., PoFa ~ 0.725. Again, this result is 
somewhat strange, since the logarithmic power 7 was left 
adjustable in the VI2 model. 

Completely analogous results were found when com- 
paring the V4 and V5 models from Table 3 . We obtain 
Pv4i = 0.616 for fixed 7 = 2 and Pysi = 0.595 when 
FMS let 7 float for the low energy data — again, it is very 
strange that letting 7 be fit by the data gave a lower 
probability than fixing 7 to be 2. 

To illustrate this anomaly, we recall to the reader that 
the difference between the V4i and the Direct Fit (DFi) 
models was that in the V4 model, the Regge powers 61 
and &2 were fixed at 1/2, whereas they were allowed to 
vary in the Direct Fit model, raising the probability from 
Pv4 — 0.616 to Pv6 = 0.721, as expected — the exact 
opposite of the V4 to V5 effect, where the logarithmic 
power 7 was varied from 2. Clearly, we question their 
minimization program, or their use of it. 



Thus, we conclude that there is no statistical evi- 
dence given in FMS that supports the conclusion that 
7 > 2, and thus, no basis for concluding that the Frois- 
sart bound is violated. 

Prediction of the 7 TeV total pp cross section — Since 
Table 1 of FMS does not contain the Totem point, its 
data could have been used to make predictions of the to- 
tal cross section at 7 TeV to compare with the Totem 
value, since it uses only lower energy data, although 
the FMS paper has not explicitly made any predictions. 
However, their plots of their results in Tables 1 and 3 al- 
low us to do so, visually. Inspection of the curve for the 
V4 model in Fig. 5 of FMS (the dash-dotted curve) shows 
that the central value of the 7 TeV prediction (from the 
data using only the lower energies, labeled = 1.8 

TeV) goes slightly inside the lower error bar of the plot- 
ted Totem result. To obtain their 7 TeV cross section 
prediction, together with its error, from their V4 model 
of Table 3 using their parameters in a standard error eval- 
uation, we find (Jtot = 96.2 ±4.5 mb, which is in excellent 
agreement with the experimental value of 98.1 ±2. 3 found 
by Totem. We see indeed that the Totem value is well 
predicted using 7 = 2, with the FMS numbers that they 
themselves calculated. 

Similar conclusions may be drawn from the Direct Fit 
model of Table 1. Using the Direct Fit model parameters 
of FSM's Table 1, we numerically calculate that the DF 
model predicts 95.4 ± 3.7 mb when only the la diagonal 
error due to the coefficient /3 is taken into account ; it 
becomes 95.4 ± 8.8 when their diagonal error due to Sh 
is included), so that their DF model prediction is also in 
good agreement with the Totem cross section, 98.1 ± 2.3 
mb. 

As in the preceding Section, we find that the FMS 
7 = 2 fits at low energy, after allowing for errors in the 
predictions due to the statistical errors in the fitting pa- 
rameters, successfully predict the Totem total cross sec- 
tion at 7 TeV, thus negating the necessity for considering 
a violation of the Froissart bound. In simpler words, the 
FMS fits are consistent with a saturated Froissart bound 
when the Totem point is included. 

Additional experimental evidence for a saturated Frois- 
sart bound — Block and Halzen Q (BH) have recently 
shown that the proton asymptotically becomes a black 
disk as s — > 00, using analyticity constraints that an- 
chor their fit at low energy to predict ultra-high energy 
experimental data at the LHC at 7 TeV, as well as the 
Pierre Auger Observatory results at 57 TeV. The BH 
model parameterizes the even and odd (under crossing) 
total cross sections and p-values and fits 4 experimental 
quantities, the 2 total cross sections (Jpp{v),app{v) and 
the 2 p-values, Pppiv) and ppp{i>), to the high energy an- 
alytic (complex) amplitude parameterizations (3], which 
guarantee analyticity and are much simpler to compute 
than derivative dispersion relations. Using v as the lab- 
oratory energy and denoting m as the proton mass, they 
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find that 



mJ 



Co + Ci In 



+C2lnM- , 
Vm/ 

Vm/ 



(4) 
(5) 



X l^ci +C27rln(— ) 



-/3,,cot(^)(ii) 
2 vm/ 

^ 2 ^ Vm/ 



47r 



/+(0) 



(6) 



using the upper sign for pp and the lower sign for pp; 
i5, a, /^p', /X, Co, ci, C2 and /+(0) are real constants, with 
/+ (0) being the singly subtracted dispersion relation con- 
stant. (T°(i/), the even (under crossing) amplitude cross 
section of Eq. Q effectively becomes 



Co + ci In 



2m2 



C2 In 



2m2 



(7) 



so that asymptotically, cr'^ — > In^ s. Thus, the high energy 
functional form used by FMS, Eq. ([2]), when 7 = 2, is 
completely equivalent to the functional form of the above 
Eq. ^ . In other word, BH replace the 3 FMS parameters 
a, j3 and Sh by the 3 linear parameters co, Ci and C2. 
Thus, the BH fit uses approximately the same physical 
assumptions as the V4 model in Table 3 of FMS, fitting 
data for cms energies y/s < 1800 GeV. 

The results [10[ of a global fit to both pp and pp 
p-values and total cross sections in the energy range 
6 < -/s < 1800 GeV, with 187 datum points and only 
5 free parameters, (5, a,ci,C2 and /+(0), is shown in Fig. 
[TJ As seen from Eq. ([6]), p — ^ as s — ?> 00, which is a re- 




FIG. 1: Proissart-bounded analytic amplitude fits to p, the 
ratio of the real to the imaginary portion of the forward scat- 
tering amplitude, vs. ^/s, the cms energy in GeV, taken from 
BH [lO[. The pp data used in the fit are the (red) circles and 
the pp data are the (blue) squares. 



quirement for a black disk at infinity. However, the tiny 
change in p from 0.135 at 1800 GeV to 0.132 at 14000 
GeV implies that we are nowhere near asymptopia, where 
p = 0. 

The fits for the pp and pp total cross sections, that 
use only 4 of the same parameters that were used for 
the p-value fit of Fig. [1] are shown in Fig. [2] The 
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FIG. 2: Froissart-bounded analytic amplitude fits to the 
total cross section, atot, in mb, for pp (dashed curve) and pp 
(dot-dashed curve) from Eq. ([5]), vs. ^/s, the cms energy in 
GeV, taken from BH . The pp data used in the fit are the 
(red) circles and the pp data are the (blue) squares. The fitted 
data were anchored by values of af^^ and crf^^, together with 
the energy derivatives da^^^/dv and da^^^/dv at 6 GeV using 
FESR, as described in Ref. We note that their ultra- 

high energy total cross section predictions that are made from 
their analytic amplitude fit use only total cross section data 
that are in the lower energy range 6 < ^/s < 1800 GeV. 

dominant In^ s term in the total cross section (see 
Eq. (O) saturates the Froissart bound [ll; thus it controls 
the ultra-high energy behavior of the total cross sections. 

A very important role in fixing the overall fit is played 
by Finite Energy Sum Rules (FESR) [ll|, ll2| anchoring 
the fit at the low energy end, taken to be y/s = 6 GeV. 
The FESR allow BH M 



to fix the fit to both crf^^ and 

pp 



crfot , together with their two energy derivatives da^^^ / dv 
and da^^^/dv, at the low end, y/s = 6 GeV, by using the 
many precise low energy total cross section measurements 
between ^/s of 4 and 6 GeV. These FESR constraints 
then can be used to evaluate directly the values of cq 
and /?-p' , two of the four parameters needed to determine 
cr°, the even high energy total cross section of Eq. (jj]). 
These fixed values of co and /3-p' , together with with the 
2 globally fitted values of ci and C2 required for (T°(z/) 
(obtained from fitting simultaneously the high energy to- 
tal cross section and p measurements in the energy region 
6 < Vs < 1800 GeV), are listed in TableHl We remind the 
reader that only data in the energy region 6 < ^/s < 1800 
GeV are used in this global fit (together with the prolific 
and accurate 4 to 6 GeV total cross section data used 
for the 6 GeV low energy 'anchor points'). We note that 
C2, the coefficient of In^(s), is well-determined, having a 
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statistical accuracy of ~ 2%. Thus, we see from Fig. [2] 
that the experimental data show that a saturated Frois- 
sart bound model is accurately satisfied for CTtot , the total 
cross sections for both pp and for pp in the energy interval 
6 < y/s < 1800 GeV; this accuracy of prediction mainly 
results from the use of the FESR constraints on the high 
energy analytic amplitude fit 11 , 12] ■ Ultra-high energy 
total cross sections, for which there are no distinction be- 
tween pp and pp interactions — both being given by cr" — 
are now accurately predicted. For example, we obtain 
values for the total pp cross section of fi" = 95.4± 1.1 mb 
at 7 TteV [13 and 134.8 ± 1.5 mb at 57 TeV '1% where 
the ±lcr errors are calculated from the (correlated) er- 
rors of their fit parameters Ci and C2. The (black) upper 
curves in Fig. [3] are plots of these BH [1] predictions for 
the total cross section cr" vs. -^/s, with the solid curve 
being the central value and the dashed curves the ztltr 
error curves. 



TABLE I: Values of the parameters, in mb, needed for the 
even amplitude total cross section, cr°(i/) of Eq. (j4}, taken 
from Ref. [T^j. 



co=37.32 mb. 


/3-p' =37.10 mb 


ci=-1.440± 0.070 mb. 


C2=0.2817± 0.0064 mb. 



The inelastic cross section, crPj^j, is determined by nu- 
merically multiplying the ratio of the inelastic to total 
cross section with the fitted total cross section (t°. The 
ratio of inelastic to total cross section was determined 
from an eikonal model, called the 'Aspen' model; for de- 
tails see Ref. and Ref. 0. 

Block and Halzen [idl found that 



-0.5 



aLi(i^) = 62.59 (-) 

-h0.14331n^ (—) mb 



-h 24.09 + 0.1604 Inf— ) 
Vm/ 



(8) 



valid in the energy domain, ^/s > 100 GeV. 

The lower (red) plots of Fig. [3] are the BH Q pre- 
dictions (Eq. ^ for high energy inelastic cross sections 
(Tinei, as a function of -^s. The solid (red) plot is the 
central value and the error bands corresponding to ztla 
are the lower dashed (red) curves. All of the existing in- 
elastic cross section measurements for pp, as well as the 
eight new ultra-high pp measurements, are shown. The 
agreement with experiment is excellent over the entire 
energy scale. In particular, the agreement with the new 
highest energy (57 TeV) experimental measurements of 
both CTtot and CTinei is striking. Since none of the experi- 
mental datum points in Fig. [3]were used in making these 
predictions, it is clear that the In^ s predictions for CTinci 
and CTtot are strongly supported by the existing ultra-high 
energy measurements. 




FIG. 3: Predictions for crtot and (Tinoi, in mb, vs. y/s, in GeV, 
for pp and pp, taken from Block and Halzen S\. For dtot, they 
have compared their predictions with recent pp TOTEM data 
at 7 TeV and Auger data at 57 TeV, while for a-moi, they have 
compared their results with 2.76 TeV pp data from ALICE, 7 
TeV pp data from ALICE, ATLAS, CMS and Totem, as well 
as with the 57 TeV pp inelastic cross section. The upper solid 
(black) curve is the central-value prediction for (Jtot and the 
lower solid (red) curve is the central- value prediction for tJinci • 
The dotted curves are the errors (ilcr) in their predictions, 
due to the correlated errors of the fitting parameters. 



Finally, BH [l^ determined the ratio of CTinoi(s) / ftot (s) 
as s — oo, given by the ratio of the In^ s coefhcients in 



CTjJjpi and CT , respectively, i.e. 



CTtot 



^2 
C2 



0.1433 
0.2817 



0.509 ±0.021, ass->cx.,(9) 



that is well within error of the expected value of ^ that 
is appropriate for a black disk at infinity. 
Conclusions — We find that: 



1. 



2. 



there is no sound statistical evidence put forth 
by FMS d for their conclusion that the Froissart 
bound is exceeded. Their models with In^ s that 
only use lower energy cross sections actually predict 
the Totem total pp cross section reasonably well. 

both the measured total pp cross sections and the 
inelastic cross sections are fit up to y/s — 57 TeV by 



a saturated Froissart-bounded In 
associated with a black disk. 



! behavior that is 



3. the forward scattering amplitude is pure imaginary 
as s — >■ oo, as is required for a black disk. 



4. the ratio of CTinci/ctot 0.509±0.021, as s 
compatible with the black disk value of 0.5. 



oo. 



Thus, we conclude that existing experimental evidence 
strongly supports the conclusion that the proton becomes 
a black disk at infinity, whose total cross section goes as 
In^ s as s — >■ oo. 
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